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ABSTRACT 

A quasi-three-dimensional rotor viscous code is used to 
predict high subsonic flow through an annular cascade of turbine 
blades. The well known Baldwin-Lomax turbulence model is used in 
the program. An attempt was made to implement a new turbulence 
model, based on renormalization group theory in the program. This 
was done to improve the prediction of the boundary layer transition 
on the blade surfaces .and subsequent wake development. The 
comparison of these two turbulence models with experimental data 
are presented. Pressure, velocity ratio, flow angle distributions 
and downstream wake predictions were studied using results from 
RVCQ3D (Rotor Viscous Code Quasi-Three-Dimensional) code. The 
computed results showed good agreement with experiment when 
comparing the blade surface local static pressure to inlet total 
pressure ratio at the midspan position of the annular turbine 
cascade. The computational approach used to implement’ the 


turbulence model is also described. 
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I. INTRODUCTION 


A. BACKGROUND 

The term renormalization belongs to quantum physics. It 
arises when the discrete formulation for particles is extended 
to include the case of continuous fields. Initially this is 
done by formulating the discrete problem on a lattice. In 
turbulence theory, when the continuum limit is taken, and the 
Kolmogorov spectrum hypotheses are applied in the wavenumber 
Space, the infrared divergence of turbulence theory arises. 
Renormalization is a standard method of removing the resulting 
Singularities and this procedure is described in McComb 
(Ref. 1). 

In an earlier paper [{Ref. 2], Yakhot and Orszag 
have proposed the application of renormalization group theory 
to turbulence modelling. In the RNG approach, the Reynolds 
stresses are accounted for by repeatedly recasting the 
equation of motion with the smallest spatial scales 
represented by effective larger scale motions. Following each 
Seal ization, the spectrum 1s narrowed by removal of a band 
of the highest remaining frequencies. When the process is 
complete, a set of averaged equations are obtained along with 


equations for an effective turbulent viscosity. 


Lund (Ref. 3] used an algebraic mixing-length 
model to calculate provisional values of eddy viscosity at 
each point normal to the wall. These provisional values are 
then processed by the RNG equation and the actual values of 
eddy viscosity to be used in the mean-flow calculation are 
determined. The RNG equation is such that if the provisional 
value of the eddy viscosity is below a certain threshold 
value, which 1s determined by the Heaviside function, the eddy 
viscosity will be zero. As the provisional values of eddy 
viscosity increase above the threshold, the RNG equation 
returns non-zero values of eddy viscosity that are always less 
than the provisional values. Since the provisional values of 
eddy viscosity increase with Reynolds number, the RNG eddy 
viscosity will be zero for low Reynolds numbers, increase for 
moderate Reynolds numbers, and finally asymptote to a fully 
turbulent value for large Reynolds numbers. This type of 
behavior should make the RNG approach well suited for 
transition calculations. The heart of the algebraic RNG model 
is a ramp function that switches on when the provisional value 
of eddy viscosity is above the threshold. If the ramp function 
is on, the RNG equation takes the form of a cubic equation to 
be solved for the eddy viscosity. Lund also mentioned a 
transformation of the cubic equation to a quartic equation and 
that this form of the RNG model would possibly be more robust. 

Kirtley (Ref. 4] used a similar quartic equation 


to study complex three-dimensional turbomachinery flows. As 


the RNG equation was expressed as a function of mean strain 
rate, the model was then referred to as "energy-based". In the 
procedure of solving the eddy viscosity, the shear stress was 
introduced into the quartic equation and a modified mixing 
length was applied in the wake region. Newton’s method was 
used to find the roots of the RNG equation, with the shear 


stress determined from a previous value of eddy viscosity. 


B. PURPOSE 

The aim of this thesis is to improve the viscous flow 
analysis in turbomachinery blade rows. The Navier-Stokes 
equations need to be solved and these are simplified by 
Reynolds averaging which gives rise to the Reynolds stress 
tensor in the momentum equations. Turbulence models are needed 
to account for the Reynolds stress tensor, and these represent 
the transport of momentum due to turbulent fluctuations. Good 
turbulence models have been developed for relatively narrow 
ranges of flows with excellent accuracy, but they can not be 
applied with confidence to complex flows. 

For computational efficiency and code simplicity, the 
Baldwin-Lomax algebraic turbulence model works very well for 
many two-dimensional flows. This model calculates’ the 
viscosity induced by the vorticity. However, it contains many 
empirically determined constants, which are determined from 


Simple two-dimensional shear-flow experiments. These constants 


are the near-wall damping of turbulence, transition location, 
and intermittency. In the analysis of turbine flow fields, the 
transition region of the boundary layer occupies most of the 
guction-side surface. The predictions of skin friction and 
heat transfer are important in turbine blades, therefore, less 
reliance on empiricism is desired. 

The RNG formalism, cannot as yet account for free-stream- 
turbulence intensity effects. The elimination of 
infinitesimally thin high-wavenumber bands is repeated over 
the renormalized spectrum, giving rise to recurrence relations 
which link the low-Reynolds-number regions of the turbulence 
with the high. Some constants appear in these relations but 
they are determined analytically from assumptions about the 
shape of the energy spectrum. The eddy. viscosity is a function 
of the mean strain rate. 

The object of this study is to apply the RNG turbulence 
model into the Rotor Viscous Code Quasi-3-D (RVCQ3D) [Ref .5 
and Ref.6] and the comparison of the results to those 
computed with the Baldwin-Lomax turbulence model. Quasi-three- 
dimensional flow is computed by numerically solving the thin- 
layer Navier-Stokes equations using a multi-stage Runge-Kutta 
scheme with a spatially varying time step and implicit 
residual smoothing. 

This study entailed the computation of two cascade test 
cases. The first test case is an annular cascade of turbine 


Stator blades which was experimentally measured by Goldman and 


Seasholtz [Ref.7]) with a one-component laser Doppler 
velocimeter. They mapped out the flow at various radial and 
axial stations by performing pitchwise blade-to-blade and wake 
surveys. RVCQ3D was used to compute the midspan location and 
comparisons were made with the experimental measurements. 
The second test case 18 a subsonic cascade of turbine 
rotor blades, experimentally investigated by Hoheisel 
[Ref.8]. The wake traverse measurements were performed 
with a wedge type probe, and the boundary layer measurements 
were Carried out on the blade suction surface with a flattened 
Pitot probe. This test case was only computed with the 
Baldwin-Lomax turbulence model and is presented as an 


Appendix. 


II. FORMULATION 


A. GOVERNING EQUATIONS 

The mass-averaged, compressible form of the Navier-Stokes 
equations theoretically model turbulent flow. The equations 
are solved numerically on a body-conforming grid system, so 
these equations are transformed from the Cartesian system 
(x,y) to a body-fitted coordinate system (£,7). The thin-layer 
approximation which neglects the viscous-diffusion terms 
parallel to the body surface is applied after the 
transformation. 

Apuby Sockol, and . Prahl (Ref .9] defined the 


nondimensionalized reference quantities as follows: 








C ie ’ x jay y Be 
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The nondimensionalized, conservative form of the two- 


dimensional, time-dependent, mass-averaged, thin-layer Navier- 


Stokes equations, from Chima (Ref. 5], are as_ follows 
(omitting the overbar): 
6) 0 s 6) z se 6, : 
— (J'F,) +J— (J 1F,) +I — (J 1F,) -—J—(J1F,) =F 2 
a5 | 9) nee y) a 2) Re Sa a) 3 (2) 


where J is the transformation Jacobian given by 


Ox Ox 
gla (3) 


Equation (2) can be rewritten as follows: 
SO en eee ah ice (4) 


to which the numerical scheme is applied. 


re) pU pV 0 
po |P4Ul pp (Pure, PL , ]puven, Pp) . _|0 
0 OV 1 pvUr+é, D 2 pvV+n, D 3 0 
e (e+p) U (e+p) V 0 
0 (5) 
Mm, U,+M,V, 
FP = Ms U, tin; Vy, 
Mm, ets (u?+v?)],+(mut+mv) u,+ (mutmyv) Vv, 


where 
- A 29 2 
my, = lean: (U+p.) 
= a 
M, = SNMy(btHe) 


, (6) 
my = (nxt ny) (B+Be) 


Ub 
m, = (Nx*ny) (Fee? 
and 
e=p (CT S (ut+v?) | (7) 
ja) MG Pel) =p (u? + v?)] (8) 


(9) 


B. TURBULENCE MODEL 
1. BALDWIN-LOMAX MODEL 
This model is a two-layer algebraic eddy viscosity 
model proposed by Baldwin and Lomax [Ref.10], and has 
been extensively used for calculating turbulent flows. The 


description of the model follows. 


The eddy viscosity is calculated using 


we (ee) inner if YSV crossover (10) 
c (ph.) aun nee Y? Verossover 


where Yu.ouover LS Given by the minimum value of y for which (4)... 


= ( U,) outer ° 


The inner and outer values are given by 


Ad 


enna = piky(1-e 4°))? fal (11) 


(p,) outer i CourerP F'n KLEB (12) 


where 


yt = Pw wy (13) 


Hy 


pememeoc, Co. = 0.0269, and Fy. iS given by 


Eueee= ME le 7, 0025 Yass Vore / Posy) (14) 


The quantities y,,, and F,,, are found from 


2 | 15 
Fly) = ylal(1-e 4) oe 


where y,,, 1S the y value for which F(y) is a maximum. 


Fries LS Given by 


Foes = (145 .5( 52%) 6] 


Max 


(16) 


The term Vp-p iS given by 
Vorr ~ (/'Vq" * Vina o Wiemann een (17) 


Except for wakes, the second term is zero. 


A mixing length is introduced in the inner viscosity 


equation, as follows; 
oo (18 
l =xy(1-e 4) 


and the w term is obtained from the vorticity, as: 


Ke 4 (vorticity) : 5 (2-2) 


—_-— 19 
pe Eee (19) 
Since the coordinates are transformed from (x,y) to 
(£,n), the vorticity in the equation above becomes 
Ov ou 1 Ov Ov Ou Ou 
=—-=— == Vp) XK eet XX, 20 


From this vorticity form, the eddy viscosity can be 


calculated by combining the inner and outer turbulent 


viscosities. 


me 


2. RNG MODEL 
In Kirtley [Ref. 4], the following relation for the 


renormalized viscosity is derived: 


wile 


v= vo (1+H(—Se A,*-C,) J 
e ; (21) 
= v,(1+H(—L?-C,)] 3 
Vo 


where H(x) is the Heaviside function, which is zero for 
negative x and equal to x for all positive x. A, is the 
wavenumber corresponding to the largest length scale in the 
inertial range. The coefficient a = 0.1186 and was derived 
(Ref. 2] from the -5/3 power law for the turbulence energy 
spectrum. C, is a function of both a and the ultraviolet 
dissipation range cut-off in the spectrum of turbulence; the 
suggested range is between 75 and 200. The dissipation rate of 


turbulence is determined from the following equations: 
€e=vs (22) 
where 


cay Chae Sl Ghoy Chae 





(23) 
Ou, 2 =) 23,/9V, du,2_ 2, du, dv\2 
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Using the same transformation techniques as in the 
discussion previously, the mean strain rate squared is 


transformed from (x,y) to (&,7) as follows: 





os? (yO! well ee ee 
S sa | Yn Bg aa Bae oe 
Ly, v_y av _y du yy du), 
ge NE Gy OE * GR’ ee 
_2 , du_y du _y av, y av), 
meer 64) EE tS “n Ge MS oat 
As discussed in reference 2° the nonlinear 


renormalized viscosity equation .can be written in a cubic 
equation form: 
Vo Viobe Ge) Va =0 (25) 

Kirtley [{Ref. 4] showed the locus of possible 
allowable roots. If the maximum value of v is chosen as the 
proper root, a jump discontinuity in the distribution of »v is 
yielded. A quartic equation transformed from equation (25) was 
used by Lund [Ref. 3], as follows: 

ee AU eIGibr) — YC) Be fo (26) 


where 


= (27) 


ne, 


and a preferable length scale proposed by Kirtley [Ref. 4] is 


the following: 


L~=C,3 tanh (7%) (28) 


Ul 





with C, = 0.0845, « = 0.372, and n is the normal distance from 
the boundary, and 6 is the thickness of the boundary layer. An 
approximation of 6, was suggested by Stock and Haase 


[Ref.11], as follows: 
Oo = 136 yaa (29) 


where y,,, 1S obtained from the Baldwin-Lomax turbulence model. 
The length scale L, in the wake region of turbomachinery flows 
is determined from a correlation derived by Raj and 


Lakshminarayana [Ref.12] as follows: 
Le =min(kxs , C,b) (30) 
where 


1 


b = 8,+¢ C# 1.35 (5 +0.02)°-5 (31) 


and 6, is the average of the pressure-side and suction-side 
boundary layer thicknesses at the trailing edge, c is the 
chord, and C,, the Gaussian constant for the wake behind a 
circular cylinder, is equal to 0.169. The equations above are 
the only empiricism in the RNG model. 

Thus, the viscosity on every grid point in the flow 


field can be calculated from the solution of the quartic 
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equation as follows. First, an initial guess of the turbulent 

eddy viscosity is set equal to zero, before the flow solver 

Starts its time marching iterative procedure. The algorithm 

for the RNG turbulence model is then as follows: 

7 1. Calculate the boundary layer thickness, mixing 
length, and mean strain rate on every grid point. 

2. Calculate the shear stress to density ratio based 

on the effective eddy viscosity at the n time step. 

3. Evaluate the Heaviside function. 

4. Apply Newton’s method to solve the quartic equation 
for the effective eddy viscosity. There are two 
posibilities: 

a. If the argument in the Heaviside function is 
positive, use that value for the Heaviside 
EUNG t Tole 

b. If the argument in the Heaviside function is 
negative, then set the Heaviside function to 
zero. 

S. Update the eddy viscosity and return to step 3, 

until the solution finally converges. 

6. The converged solution is the effective eddy 

viscosity at the n+l time step. 

7. Finally reevaluate the Heaviside function again and 
a. 1f H>0O, the n+l time step effective eddy 


viscosity is used. 
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b. if H<0, the n+l time step effective eddy 
viscosity is set equal to the laminar viscosity. 
This means the turbulent eddy viscosity equals 
zero. 

8. A turbulent eddy viscosity-to-dynamic viscosity 

ratio is then calculated for RVCQ3D. 

9. Return to main program. 

After the RNG turbulence subroutine calculates the 
turbulent eddy viscosity, the value is stored in an array and 
becomes the initial guess for the next iteration calculation. 
The calculation direction is from the leading edge toward the 
trailing edge and then on through the wake. This was done so 
that the rate of growth of the boundary layer could be 
monitored. 

The RNG model mimics the transition from laminar to 
turbulent flow within a boundary layer. As the Heaviside 
function suppresses deviations from laminar viscosity until 
the length scale and mean strain grows to values above the C, 
cutoff. The same result occurs within the turbulent-to-laminar 
transition. region. Thus the RNG model should be able to 
Capture the viscous sublayer without any near-wall damping, 
while the Baldwin-Lomax model requires the addition of a 


damping term to describe the near-wall behavior. 


as 


III. SOLUTION METHOD 


A. GRID GENERATION 

A transformation was used to transform a non-rectangular 
grid in the physical plane into a rectangular uniformly- spaced 
grid in the computational plane. 

Two-dimensional body-fitted grids were generated using the 
GRAPE code (GRids about Airfoils using Poisson’s Equation). 
This code was written by Sorenson, primarily for isolated 
airfoils [Ret . Pay A modified GRAPE code 
([Ref.14] was used for the generation of periodic C- 
type grids for turbomachinery cascades. For viscous flow, the 
grids are typically finer normal to the blade than in the 
streamwise direction. 

The locations of the intersection of grid lines with the 
inner and outer boundaries were controlled to give minimum 
shear in these regions. It is desired to have near-wall normal 
grid lines at the blade surfaces. However, because of the 
large amount of turning, normal grid lines could not be 
generated over the entire blade surface and still maintain the 
grid line periodicity on the outer boundary. The presence of 
the sheared grid lines did not appear to adversely effect the 


flow solutions. 
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B. NUMERICAL SCHEME 

The code employs an explicit, finite-difference, multi- 
Stage, Runge-Kutta algorithm. The governing equations are 
discretized using a node-centered finite-difference scheme 


ee 


(Ref. 15 and Ref. 16]: 


q, = Fi? - «Ac R( Bi?) (32) 
q, = Fi? - a. At R(q,) (33) 
q, = FS” - a,Act R(q,) (34) 
q, = Fi”) - At R(q,) (35) 
I eae fy (36) 


where a, = 1/2, @ = 1/3, a, = 1/4, and 


el j j jel _ at 
a reas itl Fy 4) i (Foy Po 
(37) 


Second- and fourth-order artificial-dissipation terms are 
added to stabilize the scheme. A spatially varying time step 
and implicit residual smoothing are used to accelerate 


convergence. 
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C. BOUNDARY CONDITIONS 

At the cascade inlet, absolute total temperature Tp), 
absolute total pressure P,, relative velocity, and relative 
flow angle are specified as constants. For subsonic flows, the 
inlet conditions are updated at each iteration by 
extrapolating the upstream-running Riemann invariant R based 


on the absolute total velocity Q to the inlet, where 


(38) 





Then, using the known total temperature T, and the isentropic 
relations, the total velocity magnitude is given by 
(y=1) Roy 4 eee rie 2 Gre Cae (39) 
oT) 

Individual velocity components are found from trigonometric 
relations by keeping the tangential velocity component fixed, 
Pressure and density are found from isentropic relations. On 
the blade surfaces, for viscous flow, the no-slip velocity 


boundary condition is used. The surface pressures are compulged 


using the normal momentum equation, namely 


mer + Neer a (ns + ne ips = -pU(n,Ug + Ny Ve) (40) 
where U = V = 0 on the surface for viscous flow. 
Surface densities are found from surface pressures anda 


specified wall temperature. Both the Baldwin-Lomax and RNG 
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turbulence models incorporated Sutherland’s viscosity law to 
calculate the dynamic viscosity. 

The exit static pressure is determined by specifying the 
exit static-to-inlet total pressure ratio. All other flow 
quantities are extrapolated at the exit plane. Periodicity is 
enforced on the "C" grid by using a pseudo-point along the 
periodic boundaries on which corresponding pitchwise variables 


are located. 
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IV. RESULTS AND DISCUSSION 
Two blade profiles that were part of the full-annular ring 
of 36 blades, with C-type gridding over the computational 
domain, are shown in Figure 1. The chord length of the actual 
vane was 55.54 mm long. The geometry of a full 250 X 60 mesh 
size (150 grid lines around the airfoil, 50 grid lines on each 
Side of the center line in the wake region) is shown in Figure 


Al in Appendix A. 
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Figure 1 The geometry of Goldman’s annular cascade 
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The experimental investigation of this cascade [Ref. 7] 
was conducted with the hub-static to inlet-total pressure 
ratio maintained at a value of 0.65. This corresponds to a 
mean-radius, ideal, exit, critical velocity ratio of 0.78. In 
Chima [Ref. 17], a boundary Cond) Eien and initial 
condition input file for the RVCQ3D code was provided. The 
computational results in the next section were calculated with 
an exit-static pressure to inlet-total pressure ratio of 
0.685. The data of streamtube variation (Appendix A.) were 


provided by Goldman [Ref. 18]. 
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A. STREAMTUBE VARIATION AND GRID RESOLUTION 
1. STREAMTUBE VARIATION 

The predicted blade surface pressure distribution, 
without streamtube variation, did not agree well with 
Goldman’s experimental data (Ref. 7]. In order to obtain a 
good baseline solution, with the Baldwin-Lomax turbulence 
model, streamtube variation based on Goldman’s [Ref. 18] 
throughflow calculations was introduced. 

The effect of streamtube variation in the convergence 


rate is shown in Figure 2. The application of a varying 


cr 
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© 
aa 
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Figure 2 Streamtube variation effects on convergence history 


Ze 


streamtube resulted in better convergence as it maintained a 
near constant convergence rate, particularly after 1000 
iterations, and finally gave an order of magnitude improvement 
on the convergence after 3000 iterations. 

Streamtube variation had a significant effect on the 
prediction of the blade surface pressure distribution as can 


be seen in Figure 3. 





Figure 3 Streamtube variation effects on pressure ratio 
distribution 


Constant streamtube thickness and radial location 
versus axial distance through the cascade resulted in a poor 
prediction of static pressure. Not only were the levels of 


Static pressure quite different, but the shape of the suction 
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Side pressure distribution was also not predicted properly. 
Finally the wake distribution predicted with and 
without streamtube variation is shown in Figure 4. Once again 
a significant improvement on not only the level of velocity in 
the wake but the shape of the wake profile is obtained with 
streamtube variation. The same streamtube variation study 
uSing a grid size of 97 X 31 is shown in Appendix B, Figures 


Bl B82) and Bae 


Circumferential Direction 





Figure 4 Streamtube variation effects on wake velocity ratio 
distribution 
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2. GRID RESOLUTION 

In order to determine the best grid resolution for the 
blade-to-blade viscous flow analysis, the influence of the 
grid size on the flow characteristics was studied. The 
solutions of three different grids, 97 X 31, 200 X 31, and 250 
X 60, were compared for the analysis of Goldman's annular 
turbine cascade. 

All of the three different sized grids were run to 
3000 iterations with the same boundary conditions and a 
Courant number of 4.0. The convergence history is shown in 


Figure 5. 
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Figure 5 Grid resolution effect on convergence history 
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Initially the convergence rate of the coarse grid is 
better than the two finer grids, however this shows signs of 
instability in the latter part of the convergence history. 
Both the fine grids (200 X 31 and 250 X 60) showed similar 
convergence histories and they both achieved approximately 
three orders of magnitude convergence after 3000 iterations. 

The effect of grid resolution on the prediction of the 
midspan blade surface pressure distribution is shown in Figure 


6. 





Figure 6 Grid resolution effect on blade surface static 
pressure to inlet total pressure ratio 


216 


———E_ 


Grid refinement from 97 X 31 to 200 X 31 shows an 
improvement in blade surface pressure distribution when 
compared with experiment, particularly on the suction side of 
the blade. Further refinement from the 200 X 31 to 250 X 60 
showed little improvement in the pressure distribution, hence 
the solution with this size grid was felt to be grid 
independent and subsequently used throughout the 
investigation. The 250 X 60 grid solution also Pave the best 
resolution at the trailing edge as the "overshoot" in P,/PO,, 
was the least for this test case. 

An investigation of the grid resolution was also 
performed on the wake profile prediction. At the 153.2 percent 
axial chord position, a comparison of experimental data with 
the prediction of local velocity to critical velocity ratio is 
shown in Figure 7. 

The coarse grid clearly produced a wake profile which 
was highly grid dependent, as the profile was not smooth at 
the wake center line. Grid refinement improved this as well as 
the prediction of the free-stream velocity. However, the 
pressure side profile shows an overshoot which is not evident 
in the experimental data. The 250 x 60 grid produced the best 
shape for the wake profile, but the deficit in the wake is 
Overpredicted and the wake diffusion or growth is 


underpredicted. The discrepancy in the free-stream velocity is 
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Figure 7 Grid resolution effect on velocity ratio at 153.2 
axial chord position 


attributed to the inability of the streamtube contraction to 
account for the three-dimensional effects in the annular 


cascade. 
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B. B-L MODEL AND RNG MODEL 

In order to test the performance of the RNG-based 
turbulence model in the RVCQ3D code, identical test cases were 
run with the Baldwin-Lomax and RNG turbulence models. Also a 
fully laminar flow calculation was included for comparison. 
Unfortunately, the RNG model did not give a converged solution 
after 3000 iterations. In Figure 8, the convergence history 
shows that the RNG turbulence model was just calculating the 


vamrnar Llow condition. 
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Figure 8 Turbulence model effect on convergence history 
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The flow velocity vector around the trailing edge for the 
Baldwin-Lomax turbulence model and RNG turbulence model are 
shown in Figures 9 and 10, respectively. The velocity vector 
plot of the Baldwin-Lomax model solution shows a near- 
symmetric double recirculation region aft of the trailing edge 
in the turbulent wake. The extent of this recirculation region 
is approximately 1.2 trailing edge thicknesses in the 
streamwise direction; however, the velocity vector plot of the 
RNG model solution was asymmetric. The suction-side boundary 
layer separated because it was still laminar near the trailing 
edge and this resulted in the large recirculation region on 
the suction side which extended aft of the trailing edge to 
approximately 2.5 trailing edge thicknesses downstream of the 
blade. Separation started from about 1/5 of the chord length 
from the trailing edge on the suction side surface. This is 


also the evident in the laminar flow test case. 
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Figure 9 Velocity vector plot of B-L model 
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Figure 10 Velocity vector plot of RNG model 
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1. BLADE SURFACE PRESSURE DISTRIBUTION 

All of the following comparisons between the B-L an 
and RNG turbulence model are made for conditions after 1000 
iterations, as this is almost at the best convergence of the 
RNG model. However the B-L results could be substantially 
improved if allowed to continue to full convergence beyond 
3000 iterations. In Figure 11, the prediction of the blade 
surface static pressure to inlet total pressure ratio is 


shown. 





Figure 11 Turbulence model comparison of blade surface static 
pressure prediction 
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Figure 8 showed that the RNG turbulence model gave a 
convergence history Similar to the laminar flow solution, and 
Figure 10 showed the laminar flow behavior of the RNG model at 
the trailing edge. But Figure 11 shows that the RNG turbulence 
model produced a different pressure distribution to the 
laminar flow case. On examination of the results over the 
whole computational domain, the turbulent eddy viscosity 
around the airfoil out to the freestream, was found to be 
equal to zero. However the eddy viscosity was not equal to 
zero in the wake region. Since the mixing length formulation 
is modified in the wake region, and is different from the 
calculation in the region around the airfoil, this would 
explain why the-RNG turbulence solution was not the same as 
the aaneimeue flow solution. And because of this, turbulent flow 
was actually calculated in the wake region. The RNG solution 
gave a better prediction of the pressure distribution, 


particularly on the suction surface of the blade. 
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In Figure 13, the flow velocity to critical velocity 
ratio in the circumferential direction at 50 percent chord 


from the pressure-side to suction-side is shown. 


Circumferential Position, Theta, Deg 





Figure 13 Turbulence model comparison of velocity-to-critical 
velocity ratio at 50 percent chord 


The Baldwin-Lomax and RNG turbulence model solutions 


at the 50 percent chord position are comparable. 
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2. MIDCHORD FLOWFIELD COMPARISON 

In Figure 12, at the 50 percent chord position, the 
blade-to-blade distribution of Cp from pressure-side to 
suction-Side is shown. Here, all three test cases show 
comparable results. Notice that the RNG turbulence model is 
slightly different from the other two test cases in the near 
Suction-Sside region. This difference is likely to be the 
effect of the recirculation on the suction-side of the blade 


near to the trailing edge. 


0.4 0.6 0.8 | 10, 
Normalized Displacement in Circumferential Direction, Theta 





Figure 12 Turbulence model comparison of Cp prediction at 50 
percent chord 
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For completeness, the predicted blade-to-blade flow 
Sao Gi stributiaon im the circumferential direction from the 
pressure-side to suction-side 1s compared Eo the 


experimentally measured data, and is shown in Figure 14. 
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Figure 14 Flow angle distribution at 50 percent chord 
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3. WAKE PREDICTION 


The wake prediction with the Baldwin-Lomax and RNG 
turbulence models and the laminar flow calculation, are shown 


in Figure 15. 
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Figure 15 Turbulence model comparison of wake prediction, 
from pressure-side to suction-side 


The Baldwin-Lomax and RNG turbulence models both give 
a good prediction of the freestream velocity magnitude, but 
the Baldwin-Lomax model gives the best wake center line 


position and the velocity deficit. 
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C. DISCUSSION 

As the RNG turbulence model case did not completely 
converge after a large number of iterations, two major 
parameters in the quartic equation which are inside the 
Heaviside function were investigated. These parameters were 
the mixing length and the shear stress. Since the RNG model 
only solved laminar flow during the initial 1000 time steps, 
this meant the Heaviside function inside the quartic equation 
was always negative. 

The constant, C,, in the Heaviside function was initially 
set equal to 200. In order to reduce the contribution of the 
negative term in the Heaviside-function calculation, the 
smallest value of C,, 75, was used. : 

Figure 16 shows the comparison of the predicted mixing 
length for both the Baldwin-Lomax and RNG models. This is 
shown for the grid line 190, which is normal to the blade 
surface, close to the trailing edge (which is at line 200). 
This comparison is done at this location, because it was felt 
that the boundary layer should be fully turbulent at this 
Station as predicted with the Baldwin-Lomax model. Transition 
was determined to take place at the 165 grid line, which is 
located at approximately 60 percent chord. Both models reach 


a Maximum of L; inside the boundary layer and these values are 
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Figure 16 Comparison of B-L and RNG mixing length 
distributions normal to the suction surface 


within 30 percent of each other. If we assume that the 
Baldwin-Lomax model-derived mixing length is correct then the 
RNG model mixing length, which is based on the boundary layer 


thickness, 6, seems to be correct. 
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In Figure 17, the shear stress distribution versus normal 
distance to boundary-layer thickness ratio is shown. Here an 
obvious discrepancy in the level of prediction of 7 is 
evident, which most probably resulted in the Heaviside 
function not going positive and therefore not initiating 


turbulent flow calculations. 
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Figure 17 Comparison of B-L and RNG) shear stress 
distributions normal to the suction surface 
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An extra check can be made on the derivation of the mean 
Strain rate, equation (24). The vorticity calculated by the 
Baldwin-Lomax turbulence model at the wall, and the mean 
Strain rate at the wall around the blade surface, are shown 
compared in Figure 18. They are in agreement everywhere except 


at a few points close to both leading edge and trailing edge. 
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Figure 18 Shear strain and vorticity distribution on the 
blade surface 


42 


The same quantities plotted along the 190 grid line is 
Shown in Figure 19. For the near-wall region, the first 15 


grid points off the blade surface are identical. 
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Figure 19 Shear strain rate and vorticity distribution normal 
to the suction surface 


Since the values of computed vorticity and mean strain 
rate are exactly the same at the blade surface, this verifies 
that the derivation and coding of equation (24) are correct. 
As the vorticity at the wall is equal to the normal velocity 
gradient and the mean strain rate is also equal to this 


quantity. 
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V. CONCLUSION 

The annular turbine cascade flowfield calculation with the 
RVCQ3D code, which is based on the Baldwin-Lomax turbulence 
model, is grid dependent. Increasing the number of grid points 
around the body is a more efficient way to improve ‘the 
solution than increasing the number of points normal to the 
body. A quasi-three-dimensional calculation which includes the 
Streamtube variation information is strongly recommended, 
particular for this annular cascade test case. 

Initial indications were that the RNG model showed similar 
convergence (a logical first check) to the Baldwin-Lomax 
turbulence model for the first 1000 iteration, after which the 
RNG model diverged. The divergence of the RNG model could be 
expected as the model is highly non-linear, and most probably 
less robust than the Baldwin-Lomax model. The comparison of 
the computed flowfield at 1000 iteration, showed that the 
Baldwin-Lomax model was overall better than the RNG model. 
However, upon further investigation it was determined that the 
RNG model only computed laminar flow over the blade surface 
and turbulent flow in the wake. As the RNG model did quite 
well in the wake region, it should be further investigated as 
a possible wake modification to the Baldwin-Lomax model. 

The RNG model calculation was initiated by assuming 


laminar flow initially as described by Kirtley [Ref. 4], and 
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this resulted in the Heaviside function always being negative. 
Another method, which should be tried is to begin from the 
high-Reynolds number limit as proposed by Lund [Ref. 3] and 
use the RNG model to update the provisional eddy viscosity; 


i.e., begin with a positive Heaviside function. 
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APPENDIX A DATA FOR GOLDMAN’S CASCADE 
The complete 250 X 60 C-type grid for Goldman’s cascade is 
shown in Figure Al. Data for the_variation of the streamtube 


thickness through Goldman’s cascade are given in Figure A2. 
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Figure Al Geometry of C-type 250 X 60 grid 
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APPENDIX B RESULTS USING A 97 X 31 GRID 
Results are given for calculations of the flow through 
Goldman’s cascade using a 97 X 31 C-type grid. The effect of 
including the effects of streamtube contraction on the 
convergence history is shown in Figure Bl, the pressure ratio 
distribution around the blade is shown in Figure B2, and on 
the calculated velocity in the blade wake is shown in Figure 


Bi. 
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Figure Bl Effect of streamtube variation on convergence 
history for the 37 x ot, guid 


5G 





Figure B2 Effect of streamtube variation on static pressure 
to inlet total pressure ratio 
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Figure B3 Effect of streamtube variation on wake velocity to 
Chitved. Velocity raclo 
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APPENDIX C T106 (Ref. 19] TEST CASE 

As another possible test case for the RNG model, a cascade 
of turbine rotor blades was considered. This section deals 
with the attempt to establish a baseline solution uSing the 
Baldwin-Lomax turbulence model of the extensive experimental 
data set as presented by Hoheisel [Ref. 19]. 

The T106 test case =e subsonic turbine cascade which 
consisted of 7 blades with a chord length of c = 100 mm. The 
blades had an aspect ratio of 3. A set of experimental data 
for a Reynolds number of 5 x 10° and 0.8 percent inlet 
turbulence intensity was used in this baseline calculation. 
The geometry of the turbine blade is shown in Figure Cl, with 
a C-type grid mesh of 200 X 30. 

The RVCQ3D code and an updated version of RVCQ3D, with 
consistent non-dimensionalization throughout (by Tweedt), were 
both used to compute the whole flow field. Baldwin-Lomax 
turbulence modelling was used in each code. The boundary 
conditions, ambient pressure, total ‘temperature, pressure 
ratio, the leading edge flow velocity, and the inlet and exit 
flow angles were provided in the experimental data set. 
Unfortunately, neither RVCQ3D nor Tweedt’s code could satisfy 


the experimentally measured inlet and exit flow conditions. 


oe. 





Figure Cl Geometry of the subsonic turbine cascade, T106 


The inlet velocity was varied over a wide range of parameters 
in the experiment. However, the RVCQ3D code could not 
reproduce complete the isentropic velocity at the inlet 


boundary. 
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1. Cp COMPARISON 
In Figure C2, both codes are seen to generate a better 
distribution of Cp on pressure side than on the suction side. 
Overall, Tweedt’s code seems to give better results than the 


—a 


RVCQ3D code. 





Figure C2 Initial prediction of the blade surface static 
pressure coefficient distribution around the blade 


os) 


ad 


2. SKIN FRICTION COMPARISON 
The skin friction coefficient (C,) distribution on the 
Suction side of the blade is shown in Figure C3. C, was 
calculated from the ratio of the shear stress to the local 
dynamic pressure at the wall. Neither prediction approaches 


the experimental data. 





Figure C3 Initial prediction of the skin friction coefficient 
distribution on the blade surface 
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